Electric field in hard superconductors with arbitrary cross section 

and general critical current law 

A. Badia-Majos^'Q and C. Lopez^ 

^ Departamento de Fisica de la Materia Condensada-I.C.M.A., 
C.P.S.U.Z., Maria de Luna 1, E-50018 Zaragoza, Spain 
^ Departamento de Matemdticas, 
Universidad de Alcald de Henares, E-28871 Alcald de Henares, Spain 

(Dated: February 2, 2008) 

Abstract 

The induced electric field E{x) during magnetic flux entry in superconductors with arbitrary 
cross section $7 and general critical current law, has been evaluated by integration along the vortex 
penetration paths. Nontrivial flux motion streamlines are obtained from a variational statement 
of the critical state, which takes the form of an optimization problem on the finite element dis- 
cretization of rj. The generality of the theory allows to deal with physical conditions not considered 
before. In particular, it is shown that the boundary condition to be used for determining E is the 
knowledge of the locus E = {) within the sample. This is emphasized for anisotropic materials in 
which the electric field is not parallel to the surface. Both numerical and analytical evaluations 
are presented for homogeneous materials with different geometries: convex and concave contours, 
samples with holes, variable curvature contours, and for anisotropic samples. In the isotropic case, 
discontinuities in the electric current paths are shown to be related to changing curvature of the 
sample's surface. Anisotropic samples display the same kind of discontinuities, even for constant 
surface curvature. 

PACS numbers: 74.25. Sv, 74.25.Ha, 41.20.Gz, 02.30.Xx 
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I. INTRODUCTION 



Type II superconductors in the Critical State (CS) hold the key property of allowing 
non dissipative volume electric currents. These can be directly fed through contact pads, 
or induced by a penetrating external magnetic field. Bean's one dimensional model for 
the CS-^ is a very simple and successful tool for determining the metastable states of a 
magnetized sample. Such states are described by the spatial distribution of current densities 
and magnetic field. The operational advantage of this model is its sharp simplification 
against the use of the complete Maxwell equations set, plus an accompanying constitutive 
law E{J) for the material. However, the calculation of transitient electric fields during the 
magnetization process has attracted much interest. In particular, the knowledge of E allows 
to evaluate the local power dissipation J ■ E, which is a basic quantity if one considers the 
thermal stability of the superconducting state. 

Although E is not an explicit variable (it is customary bypassed) in the critical state 
theory of hard superconductors,— a quasistationary model for the hysteresis losses allows 
to estimate it for such case, i.e.: E ~ 6^/6t/6i. Hence, E is locally given by the amount 
of magnetic flux which has crossed the unit length per unit time. This idea was already 
used by Swartz and Bean^ in order to produce a pioneering model for magnetothermal 
instabilities. More recently, it was brilliantly applied for predicting the observed flux jumps 
in the magnetization of high superconductors.^ However, as a shortcoming of the method, 
6i must be taken perpendicular to the magnetic flux paths, which are not obvious, unless 
for very simple geometries (infinite slabs and circular cylinders). In addition, one should 
consider the general case, in which J and E are not parallel (anisotropic material), and the 
fully vectorial information is required. 

Remarkably, as demonstrated by Brandt and co-workers^'^ for the case of superconductors 
with rectangular (and related) cross sections, the exact details of the function E{x) strongly 
depend on the sample's geometry. This was shown by analytical critical state calculations, 
as well as numerically evaluated, for the power law ~ J", and experimentally tested. On 
the other hand, the relevance of the local behavior of E{x) has been recently shown^ in thin 
film experiments. Two order of magnitude variations of E have been reported, indicating 
that averaged values may be quite meaningless. 

In this work, we propose a numerical method for the theoretical evaluation of E{x), based 
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on general critical state principles^ and well suited to problems with arbitrary geometrv- 
Since the pioneering finite element oriented modelling of superconductors proposed in Ref|8|, 
a detailed analysis of transitient electric fields has been lacking. Furthermore, within our 
theory, very general critical current restrictions may be considered. In fact, the material law 
for the superconductor is given in the form J G A, with A some bounded set around the 
origin. 

— * 

For each time step, the data for calculating E are the corresponding magnetic field 
distributions before and after the changes, as determined by our variational model. Then, 
the fiux penetration paths can be derived and, in two dimensional problems, a simple line 
integration allows to obtain the electric field at every point of the sample where the profile 
has changed {E vanishes otherwise). The method maintains the advantages of the variational 
principle. It has a simple formulation for general arbitrary geometries, it allows to consider 
anisotropic models, and it has a convenient numerical implementation on the finite element 
discretization. 

By application of the previous ideas to several cross sections, we report on the complex 
E{x) structure under various circumstances. In particular, we focus on discontinuity lines, 
originated either by corners, inhomogeneities or just by changing curvature in the case of 
smooth sample boundaries. Finally, we report on some pecularities for the case that fiux 
paths are not perpendicular to the electric current fiow. This property is outstanding for 
anisotropic material in which the current fiow is forced away from the crystal principal axes, 
as opposite to the configurations considered before.-^ 

II. THEORY 

A. Mathematical Statement 

Recall that the variational interpretation of the longitudinal critical state for a long 
superconductor with arbitrary cross-section Q reads^^ 

minimize - / |-f^n+i — -f^n^ ^ ^-L > (1) 

with i^n+i the unknown magnetic field at the time layer (n + l)6t, u = gradifn+i, and A_l a 
90° rotation of the restriction set A for the current density. In fact, if the superconductor is 
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parallel to the Z-axis, u= {—Jy, Jx)/Jc- Notice that the typical range in which equilibrium 
magnetization and surface barrier effects can be neglected is assumed, i.e.: B ~ fJ^oH. 

The solution of Eq.(PQ) may be investigated by the Optimal Control theory, as shown in 
the above mentioned articles. On following this formalism, we define a Hamiltonian density 

H{Hn+l,U,P,x) =p-U - ^\Hn+l - Hn\^ , (2) 

which contains the Lagrangian density to be minimized, as well as the associated auxiliary 
momentum p. Then, the solution of the problem (superindexed by *) verifies the canonical 
system 

OT~C 

grad H*^^ = — = u\ (3) 

with u* determined by the algebraic condition of maximality 

n{H,u\p,x) >niH,u,p,x) VmgAx. (5) 

Owing to the linear dependence of 7i on m (the Lagrangian density does not explicitely 
depend on this variable), the solution may be characterized by a maximum projection rule 

p ■ u* > p ■ u V n G A_L , (6) 

and the critical state is given by a control vector lying on the boundary of the allowed control 
set {u* G doA±). Recall that the generalized concept of the boundary doA± = {dA±} \J {0} 
must be used.— Recall also that, in the case that A_|_ is a circle this leads to the standard 
critical state expression \J\ = Jc,0. 

We want to emphasize that, from the mathematical point of view, Eq.® is valid for 
any shape of the region A. Physically, this means that our theory can host a variety of 
CS models, as it was emphasized in Ref 3 In previous work, this generality has been 
exploited for simulating the quasistatic magnetization properties of hard superconductors 
under various conditions. Here, we show that this may also be applied to the transitient 
electric fields. 

B. Physical interpretation: electric field 

Below, we show that the mathematical variable p (Lagrange multiplier of the theory) 
is physically meaningful if Eq.(jl)) is interpreted as a discretized continuity equation for the 
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magnetic field. In fact, in tlie absence of either sources or sinks one may consider the 
transport equation 

OH 

— + diY Jh = . (7) 
Then, p is simply related to the magnetic field current density by 

p= -JH6t . (8) 

Let us further exploit this equivalence, which will be the basis for obtaining the transitient 
electric field when the critical state is perturbed by the external excitation. First, notice 
that Eq.(@]) may be written as (superindices will be dropped hereafter, because the optimal 
solution is assumed) 

divp = Hn+i — Hn = gradp ■ p + p divp . (9) 

Now, define the vector field e=pxz. It is apparent that one has 

rote ■ z = —{Hn+i — Hn) = —grade ■ p — e divp . (10) 

This equation is the main result of our work, as it establishes the connection between the 
variational problem solution [Hn+i{x) and p{x)] and the physical variable e. In fact, notice 
that e is nothing but the induced electric field in the evolution H^ —>■ i/n+i if appropriate 
units are used, i.e.: Eq. ffTIH) represents a convenient discretized form of the induction law 
V X E = -B. 

We call the reader's attention that the main advantage of using Eq. pOj) relies on the fact 
that it contains the full vectorial information on the vortex trajectories (given by p), and 
this allows to derive E for nontrivial cases, within the CS framework. Previous studies were 
restricted to straight line flux paths (square and circular symmetries). 

As regards the boundary conditions to be used, the constraint that E is parallel to the 
sample's surface is not a sufficient condition for determining this vector field, as it was 
suggested in RefO On the one side, this condition is only valid for the isotropic case, in 
which E II J. Furthermore, for an arbitrary shape of the restriction region A, even the 
direction of E is unknown along the surface and has to be calculated a posteriori. Generally 
speaking, E must be determined by the specification of the locus e = p = 0. In physical 
terms, this locus is defined by two kinds of lines. First, we have e = for the free boundary 
reached by the penetrating front (i/n+i = H^) . In addition, e must vanish at the points where 
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the vector field e displays sharp bends. Macroscopically, this follows from the continuity 
condition for the tangential component of E (provided B is finite, as it should be in the 
absence of flux cutting phenomena!). 

As soon as the critical profile H^^^i{x) has been obtained, one may cast Eq. (jl(J|) in the 
form 

6e = (i^n+i — Hn — e divp) 6s , (11) 

with s the arc length measured along the streamlines of p. Notice that p (flux penetration 
paths) may be obtained by combination of Eqs.Q and (jH)) in the form p{u). Thus, an electric 
field map may be obtained just by integration with starting points at the lines e = 0. 

In order to obtain an explicit analytical form of Eq. ()ll|) . one must select a given restriction 
for the current density. For instance, within the isotropic model [| J| < Jc (i.e.: A is a circle)] 
Eqs.(jHI) (jni) lead to gradi/n+i = u = p. Insofar as in this case e = eJ / Jc, and as the 
current density streamlines are strightforwardly determined by the critical state solution, 
one is just led to solve for the scalar field e{x), and one equation suffices. We get 

6e = (ifn+i -H^-e V^H^+i) 6s , (12) 

which is a quite simple expression in terms of the critical state profiles. Other selections 
of A, however, lead to cumbersome expresions and are better treated just at the numerical 
level, keeping Eq. (fTT|l as the basis. As an example, elliptic anisotropy is treated later in this 
work. 

III. APPLICATION: ISOTROPIC CASES 

Below, we give the results of the previous procedure for isotropic samples with cross 
sections in the form of (i) a square with circular holes, (ii) a star with sharp concave corners, 
and (iii) an ellipse. These examples have been investigated numerically, by application of 
Eqs.(^ and (fT^ . The actual numerical method relies on the finite element discretization 
of our variational statement, which was described elsewhere.— Basically, Eq.(PQ) becomes 
a matrix quadratic optimization problem, when the magnetic field is expanded in terms 
of nodal mesh functions. This can be solved by a number of computational mathematics 
algorithms.— In addition, comparison with simple analytical criteria has been done when 
possible. 
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A. Square cross section with circular hole 

Let us first consider an increasing external field applied parallel along the axis of a squared 
section cylinder (edge lenght 2a, and origin of coordinates at the center) with a cylindrical 
hole (radius R and center (0, ?/c)). As one can see in Fig^ the initial fronts of penetration are 
parallel to the edges, determining successive smaller and smaller squares, x = ±Xo, y = ±yo, 
and the usual sand pile profile. The intersection of two perpendicular fronts produces a line 
of discontinuity for the current density x = ±y, and a vanishing electric field.^ Recall that 
field contours and current streamlines are the same thing for the long specimen geometry. 

As soon as the flux fronts reach the hole yo = yc + R, it is filled with a uniform value of the 
penetrating field. Afterwards, the straight fronts associated to the edges are accompanied by 
a circular front x"^ + {y — yc)"^ = emanating from the hole. The intersection of the planar 
fronts from the edges and the circular one determine parabolic lines of discontinuity for the 
vector J, given by = 4i?(?/c + R — y) (vertical parabola) and {a — y^ — R)"^ + {y — yc)"^ = 
2(a — yc + R){a — yc — R ± x) (horizontal parabolae) as it is apparent in FiglH 

The electric field can be obtained from Eq. (jl2j) . by integration with starting point at 
the hues e = 0, towards the boundary. For illustration, Fig|21 displays a number of flux 
penetration streamlines obtained with our numerical method for a given field step. In 
addition, a 3D plot of e is given, showing that one gets the expected behavior near the 
edges (compare to RefQ) as well as a considerable increase around the hole. Notice that a 
nondisplayed high flow of magnetic field is concentrated around the line x = 0, a > y > 
yc + R, which connects the hole with the nearest edge, because it is the line by which the 
hole is being refilled. Actually, for the ideal case of a perfect circular hole, E diverges as 
E ~ 6^/6t/6i and 6i 0. 

Some additional remarks, concerning the physics of flux penetration in samples with 
holes, may be done by analyzing the case of a lattice of circular holes. Fig|31 shows the 
current density contours and associated e = line structure for a 3 x 3 square lattice. It 
is apparent that, soon after the vortices have reached the outer row of holes, a physical 
boundary is established (continuous line in the figure) where the flux velocity becomes 
zero. Subsequent penetration of flux towards the sample's core only takes place across the 
intersections of this boundary and the holes, which behave as a set of point sources with 
interfering circular fronts. Further flux motion barriers are estalished at the intersections 
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between vortex trajectories emanating from the holes (straigtht dashed hnes in Fig|SJ. 



B. Tetracuspid shaped cross section 

On taking the four points (±i?, ±-R), and drawing the four quarter circumferences of 
radius R closest to the origin we obtain a tetracuspid like region with sharp spikes (see 
FigE)). Being it the cross section of a long superconducting sample, we may consider the 
problem of an increasing applied magnetic field and determine the penetrating profile as 
well as the generated transitient electric field inside the sample. FigE] shows the penetrating 
magnetic field contours, as well as the electric field corresponding to the transition from the 
displayed to the subsequent critical state. 

The main property of this example is that e is a nonlinear function of space and, as 
a consequence, the induced electric charge density during the transition is not piecewise 
constant as reported in RefO for the rectangular geometry. In fact, on using a cylindrical 
coordinate system, with the origin at the center of a circumference (f.i.: (— -R, —R)) q may 
be analytically evaluated (just one quarter of the problem needs to be considered by virtue 
of the symmetry) 

eo BR^ sin e/r^ cos^ 9 , < 9 < 9* 

, 9*<9<tt/2-9* (13) 

-eo^i^^cos^/r^sin^^ , tt /2 - 9* < 9 < tt /2 

Above, the critical angle 9* refers to the intersection between the penetrating free boundary 
and the current &en(im(jf discontinuity lines (see Fig0)). In the full penetration regime, we get 
9* = 7r/4 and a piecewise continuous behavior for the induced charge density, with jumps 
from positive to negative values at 9*. 



C. Elliptical cross section 

The long superconductor with elliptical cross-section was already considered by Campbell 
& Evetts^^ as a model system for the critical state. It was illustrated that when the flux 
front reaches the nearest center of curvature of any part of the surface, further advance in 
that direction is halted, with a resulting cusp in the subsequent field distribution. Here, we 
show that our variational statement of the critical state reproduces this behavior and allows 
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an easy computation of penetration fronts, which have a by no means trivial analytical 
evaluation, unless some approximation is used. In this sense, we should mention that, 
based on the assumption that flux penetration has the shape of an ellipse, a very accurate 
analytical solution was proposed.^* Additionally, we can give a detailed study of the arising 
electric field, as well as the e = line structure. Fig|21 shows the evolution of penetrating 
flux. On using a and h for the ellipse semiaxes, the closest curvature centers are located at 
(±a=i=6^/a, 0). This points have been marked, and one can verify that the current streamline 
structure is modified as expected when flux goes through. 

The e = structure is also peculiar to this example. As one can see in FigEl the flux 
penetration streamlines are concentrated around the closest curvature centers, display an 
increasing angle along the straight e = segment and, eventually, reach the free boundary 
line at an angle of 7r/2. 

We want to notice that this example provides the explanation of a very general fact. The 
E = Q lines related to electrical current path discontinuities will arise for any contour with 
changing curvature. Obviously, this will be always the case, unles for circular cylinders. 
Furthermore, such discontinuity lines are always confined within the sample. The case of 
corners must be considered just an idealization for which the maximum curvature centers 
have reached the boundary itself. 

IV. APPLICATION: ANISOTROPIC BEHAVIOR 

As discussed in the theoretical section of the paper, the proposed method for obtaining 
the mean electric field distribution from the variational CS theory is quite general, and 
includes arbitrary cross sections as well as anisotropic material, among other features. In 
this respect, we recall that our theory was already applied to anisotropic samples,»l^ for 
the analysis of frequency mixing phenomena in rotation experiments, but restricted to the 
infinite slab geometry. 

Here, as a final application, a circular cross section with anisotropic material law 
Jll^ ^ Jl < (A . a„ ellipse) is presented (see FigB. As it was discussed in 
R,efJl5l this law within our variational interpretation is fully equivalent to a perturbation 
theory of isotropic media, which satisfies the basic requirement x J = 0, i.e.: the J{E) 
relation is free from local loops.-^ 
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The first observation in our simulations is that the vortex penetration paths within the 
sample (along p directions) are no longer perpendicular either to the surface or to the 
electric current streamlines. In particular, as shown in the figure, this means that the 
electric field generated at the surface has both tangential and normal components, showing 
a broader generality of the method against previous studies. In more detail, if one applies 
the maximum projection rule [Eq.®] of our variational principle, flux motion lines are given 
by the vector field 

^ ^ {^^dH/dx^dH/dy) ^ i~l'Jy,J.) ^^^^ 



l^^dH/dxf + {dH/dyy ^/YJy' + J J 

This equation allows to evaluate the angles between the physical quantities involved. For 
instance, the angle between E and J is determined by 

COS a = — , == . (15) 

J^iJl + 7^2) 

In this particular example, a change of coordinates allows to transform the anisotropic law 
into an isotropic one, mapping the circular section into an elliptic section, and obtaining 
a correspondence between apparance of cuspidal points for the current flow. Thus, on 
choosing a given value for 7, we can straightforwardly compare the results of this section to 
the example in Sec. Illl CI for 6/0 = 7. This property was already exploited in RefllJ. but 
restricted to the examination of the current streamlines and flux fronts. In our case, we have 
calculated the cuspidal points for the anisotropic system, which are plotted in FigEJjust for 
validation of the numerical method. More relevant from the point of view of this work is that 
anisotropic material laws and variable surface curvature geometry have equivalent effects on 
the electric field evolution. In particular we show the appearance of the inner points where 
the = lines associated to sharp bends in the current flow arise. However, contrary to the 
isotropic case, the vortex trajectories are no longer perpendicular to the current streamlines 
and E is no longer parallel to J [see Eq. fll5|) for the actual angle between these quantities]. 



V. SUMMARY 



A numerical method is presented which allows the computation of induced local electric 
fields during magnetic flux penetration (or exit) in hard superconductors. The theory is 
based upon a variational interpretation of the critical state, from which one can derive the 
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flux penetration streamlines (vortex trajectories) in nontrivial cases. Integration along these 
lines, starting at the E — points, and ending at the sample's surface, is used for calculating 
the electric field. The generality of the theory allows to investigate the physical properties 
related to an arbitrary sample's cross section, as well as the implications of considering a 
wealth of critical state models. In fact, the critical current law is formulated as the very 
general restriction J e A for the current density. 

The application of our algorithm has allowed a quantitative study of (i) the electric field 
in the vicinity of a hole near the edge of a square sample, (ii) the simulation of nested fiux 
barriers in the magnetization of samples with a square lattice of holes, (iii) the properties 
of E for samples with sharp concave spikes, (iv) the discontinuity lines for the current 
density, as related to changes in the surface curvature, and (v) electromotive forces arising 
from vortex penetration non perpendicular to the electrical current streamlines (anisotropic 
regime) . 

The wide application range of the theory allows the expeditious incorporation of new 
physical phenomena as spatial inhomogeneities in the pinning force, or fiux cutting interac- 
tions for non parallel vortex lattices (just by selecting A properly). 
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FIGURE CAPTIONS 



FIG. 1: Current density streamlines in the critical state for a square superconductor with a circular 
hole close to the upper edge. Four stages of the initial magnetization process are depicted. The 
diagonal e = lines have been marked in the upper left picture and dashed parabolic e = lines 
around the hole for the complete set. 

FIG. 2: Flux penetration streamlines (left half) and induced electric field modulus (right half) 
in the magnetization process of a square superconductor with a circular hole close to one edge. 
Parabolic and diagonal e = lines have been marked. 

FIG. 3: Current density streamlines in the critical state of a square superconductor with a (3 x 3) 

square lattice of holes. The e = structure has been emphasized by dark lines. Continuous style 
is used for the flux penetration barrier and dashed for connections with inner boundaries. 

FIG. 4: Upper left: Current density streamlines in the critical state for a superconductor with 
tetracuspid like cross section; right: induced electric field modulus contours for increasing magnetic 
field. Lower left: flux penetration streamlines; right: 3D plot of the induced electric field modulus, 
obtained by integration along the streamlines. Only one quarter is depicted for clarity. 

FIG. 5: Flux penetration in a superconductor with elliptical cross section in the critical state. 
Upper: current density streamlines previous (left) and subsequent (right) to the contact of the flux 
front with the maximum curvature centers. Lower left: detail of the flux penetration paths around 
a maximum curvature point. This point and e = lines have been highlighted. Lower right: 3D 
plot of the transitient electric field modulus for increasing magnetic field. For clarity, one half of 
the problem is shown. 

FIG. 6: Flux penetration in a superconductor with circular cross section and anisotropic critical 
current. Upper left: current density streamlines subsequent to the contact of the flux front with the 
cuspidal points (see text). Upper right: critical current distribution. Lower left: flux penetration 
paths in the neighbourhood of the cuspidal point. Lower right: 3D plot of the transitient electric 
field modulus. The arrows indicate the direction of E on the sample surface. 
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